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ABSTRACT: We study the influence of strong external magnetic fields on gluonic and fermionic observ- 
ables in the QCD vacuum at zero and nonzero temperatures, via lattice simulations with Nf = 1 + 1 + 1 
staggered quarks of physical masses. The gluonic action density is found to undergo magnetic catalysis 
at low temperatures and inverse magnetic catalysis near and above the transition temperature, similar 
to the quark condensate. Moreover, the gluonic action develops an anisotropy: the chromo-magnetic 
field parallel to the external field is enhanced, while the chromo-electric field in this direction is sup- 
pressed. We demonstrate that the same hierarchy is obtained using the Euler-Heisenberg effective 
action. Conversely, the topological charge density correlator does not reveal a significant anisotropy 
up to magnetic fields eB ~ 1 GeV 2 . Furthermore, we show that the pressure remains isotropic even 
for nonzero magnetic fields, if it is defined through a compression of the system at fixed external 
field. In contrast, if the flux of the field is kept fixed during the compression - which is the situation 
realized in the lattice simulation - the pressure develops an anisotropy. We estimate the quark and 
gluonic contributions to this anisotropy, and relate them to the magnetization of the QCD vacuum. 
After performing electric charge renormalization, we obtain an estimate for the magnetization, which 
indicates that QCD is paramagnetic. 
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1 Introduction 

Strong magnetic fields are an interesting probe of Quantum Chromodynamics (QCD). They are rel- 
evant in at least three strongly interacting physical systems: the early universe, magnetars and non- 
central heavy ion collisions [1-3]. Moreover, external magnetic fields challenge our theoretical under- 
standing of non-perturbative phenomena in QCD. The chiral magnetic effect [4, 5], where the magnetic 
field may couple to topological features of QCD (on an event-by-event basis), has triggered a lot of 
activity in this direction. 

A rather robust effect is the enhancement of the quark condensates by the magnetic field B, 
and, thus, the generation of dynamical mass, the so-called magnetic catalysis [6]. The perturbative 
approach to this mechanism (and also to other effects) focuses on the Landau levels, i.e., on the 
eigenvalues of the Dirac equation in constant external magnetic backgrounds. The lowest Landau 
level has vanishing eigenvalue and a degeneracy proportional to the magnetic flux [7]. Therefore, the 
magnetic field increases the density of small eigenvalues and thus, via the Banks-Casher relation [8], 
the quark condensate. The enhancement of the condensate at zero temperature is also directly related 
to the positivity of the scalar QED /3-function [9]. Magnetic catalysis has been confirmed in various 
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non-perturbative approaches to QCD, as well as in ab-initio lattice simulations [10-14]. For a recent 
review see, e.g., ref. [15]. At finite temperature, however, inverse magnetic catalysis was found for 
physical quark masses [12, 13], manifested by a reduction of the condensate with growing B. This 
effect reduces the transition temperature T c of the QCD crossover [12]. Inverse magnetic catalysis may 
be attributed to the sea quark back-reaction being particularly effective near the transition [16]. 

The presence of a uniform homogeneous magnetic field breaks the isotropy of space. The Landau 
level eigenfunctions are plane waves in the direction of the magnetic field, whereas they are (localized) 
harmonic oscillator eigenfunctions perpendicular to B. Accordingly, in the interacting case, the Dirac 
eigenmodes are expected to be elongated [17] along the magnetic field, and for asymptotically large 
B, the dynamics of the system effectively reduces to 1 + 1 dimensions. The anisotropy will result in 
nonzero expectation values of certain observables that at zero field vanish due to rotational symmetry. 
A prominent example is the magnetic susceptibility of the QCD vacuum [18], the proportionality 
factor between the spin polarization {tpaxyip), the external magnetic field (directed along the positive 
z axis) B = B z = F xy and the condensate {ipip)- Its value has recently been determined using lattice 
simulations at physical light and strange quark masses, in the continuum limit [19] (for quenched 
results see [20, 21]). 

In this paper, we mainly focus on gluonic quantities. Since only the quarks are electrically charged, 
the effect of magnetic fields on the gluons may at first seem secondary. Nevertheless, we find the 
total gluonic action and the anisotropic magnitudes of the field strength components to reveal a pro- 
nounced dependence on the magnetic field. Interestingly, the gluonic action behaves quite similar to 
the condensate. At zero temperature it is increased by the magnetic field, around the transition it 
is non-monotonic, whereas above T c it undergoes inverse magnetic catalysis. Note that the gluonic 
action is in some sense tied to the quark condensate: the sum of both enters the interaction measure 
(trace anomaly). We will exploit this relation to improve the definition of the gluonic action at our 
nonzero lattice spacings. 

Our second objective is to determine the anisotropy between the (squared) gluonic field strengths 
parallel and perpendicular to the external field. This anisotropy was first measured on the lattice 
by the Berlin group [14], for SU(2). The following hierarchy was found: the chiomo-magnetic field 
parallel to B has a larger, and the chromo- electric field parallel to B a smaller expectation value than 
the components perpendicular to the magnetic field. The latter all share the same expectation value at 
zero temperature due to the remaining symmetry (at nonzero temperatures, an additional anisotropy 
is induced between the chromo-magnetic and chromo-electric fields). We confirm this hierarchy for 
real QCD with three quark flavors of physical masses. For a comparison, we give the Euler-Heisenberg 
effective action for small color fields in a constant external magnetic field background. It reveals the 
same hierarchy. 

We also determine the anisotropics of the fermionic part of the QCD action, and discuss how 
these - together with the gluonic anisotropics - are related to the diagonal components of the energy- 
momentum tensor. The subtle question, whether the hydrodynamic pressure in the presence of mag- 
netic fields is isotropic or anisotropic is the topic of an ongoing discussion, see e.g. refs. [22-26], and 
is discussed in sec. 4. This question is highly relevant, for instance, for density fluctuations in the 
early universe [27], for magnetized neutron stars [2] and for the elliptic flow observed in heavy-ion 
collisions [28]. We show that the pressure is isotropic, if defined via the compression of the system 
with the magnetic field kept fixed, whereas it is anisotropic if the flux of the field is kept constant 
during the compression. We demonstrate how the latter anisotropy can be used to obtain the mag- 
netization of the QCD vacuum, and give the first estimate on the renormalized magnetization from 
lattice simulations. 
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For the third gluonic quantity under study, the two-point correlator of the topological charge 
density, we do not find a significant anisotropy for the applied range of magnetic fields. One might 
have expected an anisotropy mediated by elongated quasi-zero modes, given the strong link between 
topology and chirality. Therefore, also this null result contains interesting information about the 
response of the QCD vacuum to a magnetic field. 

This paper is organized as follows. In the next section we define our observables in the continuum 
and on the lattice, and discuss their renormalization. In sec. 3 we briefly describe our simulation 
setup and present our main numerical results. Subsequently, in sec. 4, we address the anisotropy 
in the pressure and relate this to the QCD magnetization. This is followed by our conclusions in 
sec. 5. The article is augmented by four appendices that contain more detailed derivations and the 
Euler-Heisenberg effective action for our setting. 



2 Observables 

On the lattice, observables are given in terms of the partition function, which is written as the func- 
tional integral over the gauge links U, 

Z= I VUe- ps " JJdet [a#($) + am f ] , (2.1) 
/ 

where f = u,d,s labels the quark flavors, a is the lattice spacing, and the parameters are the inverse 
gauge coupling j3 = Q/g 2 , the lattice quark masses amf and the flux $ of the external magnetic field 
eB, which couples to the quark determinants through the electric charges qf/e of the quarks, see 
below. Here e > is the elementary charge. The quark condensates and the gluonic action are given 
by the partial derivatives with respect to the lattice parameters, 



d(am f y xr/r// NlN t x J J/ ' d(3 ' JVfM 
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(2.2) 

where we also defined the corresponding densities in physical units. Here, N s and Nt are the number of 
lattice sites in the spatial and temporal directions, respectively, and a is the lattice spacing. In terms 
of these, the three-volume and the temperature read V = (N s a) s and T = (iVja) -1 . For convenience, 
we define for any observable O, 

AO= (0)\ eB - (O)\ , (2.3) 
as the difference of expectation values with and without the external magnetic field. 

2.1 Interaction measure 

The gauge action density and the quark condensate densities enter the interaction measure (trace 
anomaly) , 

TdlogZ . . 

an observable relevant for the QCD equation of state. We consider the gluonic contribution to AJ, 
which - after the improvement procedure discussed in app. B - reads 
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The first term on the right hand side of eq. (2.5) is the direct gluonic contribution, see, e.g., ref. [29]. 
Since the additive divergences are canceled in As g (see the discussion in apps. A and B), this first 
term already has a well-defined continuum limit (at B = 0, typically differences in the temperature 
are utilized to achieve this). The leading lattice discretization error, however, is of 0(1/ log a) and, 
therefore, vanishes very slowly. The second, fermionic term in eq. (2.5) does not contribute to the 
continuum limit, but it cancels the logarithmic discretization error, thereby improving the convergence 
to C(a 2 ), see app. B. The improvement is evident from our numerical data, see fig. 7. The procedure is 
equivalent to multiplying the result by a finite renormalization constant 1 + 0(1/ j3) that we determine 
non-perturbatively, using the line of constant physics (LCP) for the action that we employ [29]. To 
our best knowledge, such an 0(l//3)-improvement of a gluonic quantity by a fermionic quantity has 
not been considered in the literature previously. 

2.2 Anisotropies 

The continuum counterpart of the gluonic action in Euclidean space is written as 

PSgt+^J d 4 xtr J^(x) = 4/ d 4 xtr [S 2 (x) + B 2 (x)} . (2.6) 

The field strength is defined in terms of the SU(3) gauge potential as = d^A u —d u Afj,+i[Afj,, A u ], 
and consists of chromo-electric E{ = and chromo-magnetic components B{ = €ijkJ~jk/^- We use 
calligraphic letters to denote the non-Abelian SU(3) fields, to distinguish these from the external 
Abelian field 1 . The full covariant derivative reads D„f = d^ + iA^ + iqfA^. Without loss of generality, 
we will take the external magnetic field B to point in the z-direction. The simplest gauge field to 
realize this (in infinite volume) is A y = Bx, A^ = (// ^ y). 

The translation of these quantities to the lattice discretization is straightforward. For the gauge 
action we use the tree-level improved Symanzik action [30], 

s 9 = S? m = E E \ Retr P K»(*)> ( 2 -7) 

where -P^(n) denotes a sum of gluonic loops lying in the [i-v plane, see eq. (C.4), and n runs over 
lattice sites. Therefore, S g is readily decomposed into planar components and, therefore, into squared 
traces of the chromo-electric and chromo-magnetic field strengths, according to eq. (2.6), 

tr£ 2 (n) = 2 Retr P 4i (n), tr£ 4 2 (n) = 2^ \e ijk \ RetvP jk (n). (2.8) 

j<k 

In the following, the components in the direction of the external field B [| z are denoted as parallel, 
whereas the x and y components as perpendicular, 

£ 2 +£ 2 B 2 + B 2 

£ 2 = £ 2 , B\ = Bl £i = ^p, B l = ^y. ( 2 . 9 ) 

We define the anisotropies as the densities of the expectation values of differences between these 
components, 



A ^ = % (f E ( tr ^ n ) " tr£ t ) • A &) = y ( f E (t^Kn) - trB 





2.10) 



1 Note that the Euclidean £ 2 (E 2 ) turns into — £ 2 {—E 2 ) in Minkowski space-time, whereas the sign of the squared 
magnetic fields remain the same. 
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which are intensive, gauge invariant quantities. The factor of /3/6 = g 2 is included here to define the 
anisotropy of the total action, eq. (2.6). At B = 0, the spatial isotropy of the system is restored, and 
both A(S) and A{B) vanish. 

Now we turn to the anisotropy induced in the fermionic part of the action, ^ flfi'fy f. Even though 
the Dirac structure of ij) in the staggered discretization is related non-trivially to the continuum 
structure, the fermionic action can still be unambiguously separated into four terms, corresponding to 
space-time directions [i. Each quark flavor / gives a contribution, 

<W = (2.H) 

where the Jj) „ are defined in eq. (C.5) for staggered quark actions. In analogy to eqs. (2.9) and (2.10), 
we define the parallel and perpendicular components and the anisotropy as 

C U = C zJ , C.j C ^[ 2 C ^. A(C f ) = T (C±J -C lf ). (2.12) 

Again, A(Ct) vanishes in the absence of the external field. In sec. 4, we will discuss the renormalization 
properties of these anisotropies and their relationship to the QCD magnetization. 

2.3 Topological charge density and its correlator 

The density of topological charge in the continuum reads 

1 1 

1\ x ) = ^^ P <jtTj r ll u{x)J r pa{x) = — y tr£(x)B(x) . (2.13) 

This definition is not suited for a direct application on the lattice, as the latter is discrete, and short 
range quantum fluctuations would dominate the signal. Therefore, we apply 5 steps of improved 
stout smearing [31] (filtering) with the standard smearing parameters p = 0.06 and e = —0.25. The 
3-loop-improved field strength tensor [32] is then measured and used in the above formula to define 
q(x). This technique has been applied in various studies of topological properties [33, 34], since it 
gives a topological density similar to that from the fermionic definition containing all eigenmodes of 
the overlap operator [35, 36]. Less filtering leads to ambiguous results for the improved field strength, 
whereas stronger filtering eventually washes out the topological structures. 

The topological density is a CP-odd variable and, thus, vanishes on average in the QCD vacuum. 
This remains true in the presence of an external magnetic field (only in combination with an additional 
external electric field E a second CP-odd variable exists, EB, that could couple to q, see, e.g., refs. [37- 
39]). 

We will consider the two-point correlator of the topological density, (q(x)q(x')) . It is a CP-even 
quantity that, at B = T = 0, depends only on the absolute value of the distance of the two points, 
due to translational and rotational invariance. At B > 0, rotational symmetry is lost, and therefore 
the correlator may be different if the four-dimensional difference vector x — x' is perpendicular to the 
magnetic field or parallel to it. At T > 0, the only parallel direction is the z-direction, whereas at 
T = 0, the parallel directions also include the ^-direction, 



{q(0)q(r)) ± 

(q(0)q(r))\\ = (q{x)q(x')} with \x — x'\ = r and x — x' 6 < 



(x, y)-plane 

(z, (t))-plane . (2.14) 
arbitrary 



We consider the difference of these correlators to quantify the anisotropy in the topological sector. 
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The integral of the full correlator over r gives the topological susceptibility T/V ■ (Q 2 ), and is thus 
positive. At intermediate distances, however, the correlator becomes negative, with the position of the 
zero proportional to the lattice spacing. This is clear from the continuum limit, in which the correlator 
is negative for all distances apart from zero, where it contains a contact term [40]. The scaling of the 
position of the zero has turned out to be very similar for many current fermion discretizations [34]. 
We will focus on intermediate distances r/a = 2 . . . 4. Note that, due to the restriction of r to 
two-dimensional hyperplanes, less discrete distances are available for the perpendicular and parallel 
correlators than for the full propagator, resulting in lower effective statistics. 

Finally, we remark that the use of smearing and of the improved field strength definition may affect 
the anisotropy of the correlator, as both techniques effectively amount to an averaging over space-time 
regions in a spherically symmetric way. For the mild averaging that we employ, however, these regions 
do not overlap strongly and, thus, the correlators defined above still contain the information about 
anisotropics in the topological charge. 



3 Results I: gluonic and fermionic observables 

Our measurements have been performed on the same configurations as used in our previous studies 
of magnetic fields in QCD. The configurations at zero and nonzero temperature have been generated 
with the tree-level improved Symanzik gauge action and Nf = 1 + 1 + 1 stout smeared staggered 
fermions, at physical quark masses, for details see [12, 13, 19]. The light quark masses are set equal, 
m u = rrid = Tn u d, whereas the strange quark mass is m s = 28.15 • m u d- The quark masses are tuned as 
a function of /3 along the line of constant physics (LCP) [29], which ensures that for all lattice spacings, 
the hadron masses are at their physical values. The quark charges are —q u /2 = q^ = q s = — e/3. 

3.1 Interaction measure 




Figure 1. The change in the gluonic contribution to the interaction measure, eq. (2.5), at zero temperature 
(left, including the continuum limit from four lattice spacings) and the corresponding light quark contribution 
m u d^-ipud'4 , ud (five lattice spacings and the continuum limit), in the same units. 



We start the analysis by considering the change of the renormalized gluonic action, i.e., of the glu- 
onic contribution to the interaction measure — A/g mp of eq. (2.5) at zero temperature. We use four 
different lattice spacings with magnetic fields eB up to about 1 GeV 2 , and perform a combined spline 
interpolation and continuum extrapolation to obtain the a — > limit. The results, together with the 
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continuum limit, are shown in the left panel of fig. 1. Clearly, — A/g mp is enhanced by the magnetic 
field, which we interpret as magnetic catalysis for gluons. For large magnetic fields, the initial quadratic 
dependence seems to turn linear, just like for the quark condensate in the same system. Note that the 
fermionic counterpart of — AI^ mp is — A/y mp = Y2f m -j '/V '} '> see the discussion in app. B. Therefore, 
for comparison, in the right panel of fig. 1, we also display the change in the average light condensates, 
m u d^'ipud' t Pud = m ud (A^Vm + A-'ipdipd) /2, taken from ref. [13], and plotted in the same units as the 
gluonic observable in the left panel. Note that the relative errors are larger in the gluonic quantity, 
since here 0(a -4 ) divergences cancel in the difference A, whereas the leading divergence for the quark 
condensates is only of 0(a~ 2 ) [29]. 

In absolute terms, the effect of the gluonic magnetic catalysis is by two orders of magnitude larger 
than that of the light quarks. Notice however, that the natural scale of the quark condensate at 
T = B = is given through the Gell-Mann-Oakes-Renner relation by M%F 2 /2 m 6 • 10~ 5 GeV 4 with 
M n being the pion mass and F the pion decay constant, whereas conventional estimates for the gluon 
condensate are of the order of 10~ 2 GeV 4 [41]. Therefore, the relative effect of the catalysis is stronger 
in the quark sector, as expected, since quarks experience the primary effect of the external magnetic 
field. We remark moreover, that the strange contribution is even larger, as it comes with the factor 
m s = 28.15 • m u d- 




Figure 2. The change in the gluonic contribution to the interaction measure at temperatures around the 
transition (left, continuum limits). The magnetic catalysis turns into inverse magnetic catalysis, very similar 
to what has been found for the quark condensate [13] (right panel). 

Next, we turn to the transition region and perform a similar combined continuum extrapolation 
of — A/g mp as we did at T = 0, this time using three lattice spacings with Nt = 6, 8 and 10. We find 
the behavior of — A/g mp to become quite different, when the temperature approaches T c ~ 150 MeV, 
see the left panel of fig. 2. Slightly below this temperature, the dependence of the gluonic contribution 
to the interaction measure on the magnetic field starts to become non-monotonic. At even higher 
temperatures, — A/g mp decreases monotonically, with magnitudes comparable to those at zero tem- 
perature: the gluons undergo inverse magnetic catalysis. Again, this behavior is very similar to that 
of the light quark condensate, as found in ref. [13] (see the right panel of the figure), with the notable 
difference that gluonic inverse catalysis persists for high temperatures, whereas for T > T c , the quark 
condensates eventually approach zero, due to chiral symmetry restoration. 
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3.2 Anisotropies 

Here, we study the individual components of the gauge action, as given by eq. (2.9). We remark that 
for the anisotropies - unlike in sec. 3.1 above - we do not perform the continuum limit, but only show 
the scaling tendency of the results with the lattice spacing. In order to carry out a proper continuum 
extrapolation, one has to subtract terms ~ (el?) 2 log a that arise from charge renormalization (see 
app. A). We will revisit this issue in sec. 4. 
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Figure 3. The components T/V (J2 n tr£>f(n)) and T /V (J2 n tT Ef (ri)) in lattice units a -4 , as measured on 
a 24 3 x 6 lattice at a temperature T = 189 MeV. The anisotropies induced by the temperature and by the 
magnetic field are indicated by the arrows. 

First, we demonstrate the hierarchy of the gluonic components at T > T c and B > 0, where effects 
from both the temperature and the magnetic field are present. In fig. 3, we plot the expectation values 
of the densities of the individual components eq. (2.9), as determined on our Nt = 6 lattices. In the 
absence of the magnetic field, the anisotropy is induced solely by the temperature, separating the 
chromo-magnetic and chromo-electric components. For B > 0, in addition the parallel and perpen- 
dicular components split, due to the spatial anisotropy induced by the magnetic field 2 . The generated 
hierarchy is (tiB"^ > (trB^J > (trSY) > (tr£,p, similar to what was observed in the SU(2) theory 
in ref. [14]. 

To determine the dependence of the anisotropies on the external magnetic field, in the left panel of 
fig. 4, we plot A{£) and A(B) (see their definition in eq. (2.10)) as functions of eB at T = 0. The parallel 
chromo-electric field is suppressed with respect to the perpendicular fields, resulting in a positive 
A(£ ), whereas the chromo-magnetic sector shows the opposite effect, giving a negative A(B). This 
non-perturbative finding is in-line with a perturbative treatment of the anisotropy, see the generalized 
Euler-Heisenberg calculation in app. D, in particular eq. (D.5). According to this calculation, trf| 
increases the effective action (to bi-quadratic order in B and J 7 ), and is thus suppressed. In contrast, 
tri3 2 reduces the action, and is favored. This implies A(£) > and A{B) < 0, as we have found. 
Furthermore, within the present statistical accuracy, the two anisotropies have the same magnitude. 

The gluonic anisotropies do not show any significant finite volume effects, and we also find these 
to be roughly independent of the temperature up to our largest T = 189 MeV. However, the signal- 

2 Note that the data in fig. 3 contains an additive divergence ~ a~ 4 . The anisotropies induced by T and B (indicated 
by the arrows in the figure) are, however, ultraviolet finite. 
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Figure 4. Anisotropics in the squared field strengths, eq. (2.10), (left panel) and in the fermionic action 
eq. (2.12) (right panel) at zero temperature. 

to-noise ratio becomes worse at high temperatures (which, in our fixed Nt approach, correspond to 
finer lattices), again due to the cancellation of 0(a -4 ) divergences in A{£) and A(B). 

The fermionic anisotropics defined in eq. (2.12) also develop nonzero expectation values for B > 0, 
see the right panel of fig. 4 for our zero temperature results for the up quark. We find the anisotropies 
A(Cf ) to be negative for all three quark flavors. The anisotropy in physical units is by about a factor of 
five larger than the anisotropies found in the gluonic sector. Similarly as for the gluonic anisotropies, 
we find the magnitude of A(Cf) to be roughly independent of the temperature. We stress again that 
the anisotropies presented here are still subject to an additive renormalization, which we discuss in 
sec. 4 below. 

3.3 Topological charge 
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Figure 5. Left panel: correlator of the topological charge density, eq. (2.14), at eB = 1.1 GeV 2 at vanishing 
temperature on a 40 3 x 48 lattice of lattice spacing a = 0.1 fm. The perpendicular and parallel (green and 
red points, respectively) correlators are compared to the total one (dashed line). Right panel: the difference 
between the total correlator at eB = 1.1 GeV 2 and at eB = (red triangles), and the anisotropy between the 
parallel and perpendicular correlators at eB = 1.1 GeV 2 (blue squares). 
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We proceed by investigating the correlator, eq. (2.14), of the topological charge density. In the 
left panel of fig. 5, the parallel, perpendicular and total correlators are shown for our largest magnetic 
field, eB = 1.1 GeV 2 on a 40 3 x 48 lattice with a = 0.1 fm, at T = 0. The results exhibit the typical 
dependence of the correlator on the distance, with a positive peak at r = 0, a zero at around r = 2a, 
and a negative tail for large distances. To quantify the effect of the magnetic field, in the right panel 
of the figure, we plot the difference between the parallel and perpendicular correlators (squares), and 
that between the B > and B = correlators (triangles). The latter reveals an enhancement of the 
correlator for small distances as B grows. Note that due to the presence of the contact term, the 
correlator at r = diverges in the continuum limit (approximately as a -6 , as observed in ref. [34]). 
For fine lattices, this divergence dominates the difference created by the magnetic field, resulting in 
a lower signal/noise ratio for r < 2a. Altogether, we conclude that the magnetic field generates an 
0(1O~ 7 GeV 8 ) upward shift in the correlator at r = 0. 

We remark moreover that there is a slight tendency of a suppression of the parallel correlator with 
respect to the perpendicular one around r = 0.2 — 0.3 fm, which is however, at the present level of our 
statistics, consistent with zero. Thus - in sharp contrast to the pronounced effect of B on the gauge 
action and on its components - we conclude that the magnetic field induces no strong anisotropy in 
the topological charge density correlator, as was also conjectured in ref. [17]. 

4 Results II: Isotropic and anisotropic pressures 

We are now in the position to discuss the isotropy properties of the pressure for nonzero magnetic 
fields. It is well known that the pressures are related to the spatial diagonal components of the 
energy-momentum tensor T a ^ [42]. In the presence of external magnetic fields, T a @ has been studied 
extensively in the last decades. In several studies, the pressures obtained from these components were 
found to be anisotropic (as, e.g., in refs. [23, 26, 43, 44]). Such a difference between the longitudinal and 
transverse pressures would have important consequences in, e.g., the physics of magnetized neutron 
stars [44]. However, in the above studies, the currents generating the magnetic field are not taken 
into account, and thus, the considered system is open, with an energy-momentum tensor that is not 
conserved [45, 46]. If the system is defined to include these currents as well, additional contributions 
arise [24, 47], which are the topic of an ongoing discussion [24, 25]. Here, we address the problem 
from a somewhat different aspect, and consider the pressures as the response of the thermodynamic 
potential of the system against compressions in the corresponding directions. 

4.1 In the continuum 

Let us consider a polarizable medium in a finite volume V = L x L y L z , at temperature T = 1/Lf, 
exposed to a magnetic field eB in the positive z direction. The magnetic flux through the x — y plane 
is given as $ = eB ■ L x L y . The partition function of the system depends on the spatial extents, the 
temperature and the magnetic field, log Z(Li, T, eB), with i = x,y,z. The energy density e and the 
pressures pi are given as total derivatives with respect to the corresponding extents, 

ldlogZ T dlogZ 

The reaction of the system to the simultaneous rescaling of all extents £„, is related to the interaction 
measure (trace anomaly) /, which we have already defined in eq. (2.4). In terms of the energy 
density and the pressures it reads I = e — p x — p y — p z . The Independence of log Z is related to the 
magnetization, 
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It is crucial now to distinguish between two distinct settings, which can be used to define the 
transverse pressures {p x and p y ) and the interact ion measure. For the derivatives with respect to L x 
and L y , one can keep fixed either the magnetic field (eB = const.), or the magnetic flux (<£ = const.). 
We will refer to these cases as the i3-scheme and the <3?-scheme, respectively. In the former case the 
flux grows as L x (or L y ) is increased, whereas in the latter case, eB decreases for increasing L x (or 
L y ). Along constant eB, the only dependence of log Z on the extents is the explicit one, whereas for 
the $-scheme, an additional term emerges due to the implicit dependence of B on L x and L y , 



(B) T dlogZ m T dlogZ TdlogZ dB 
Pi = 77 ETT— T ■ P = 77 7777777- + 



V dlogW * % V d\ogU V dB dlogU 



(4.3) 



revealing the relation between the pressures in the two schemes. The energy density - just as p z - is 
independent of the scheme, 

pf ] =P?\ P { $=P { JS-M.eB, e(*W B ), (4.4) 

which makes the interaction measure also scheme-dependent, 

j(*) = j(B) + 2M • eB. (4.5) 

Below we will omit the superscripts from the scheme-independent quantities e and p z . Note that at 
zero temperature, the z and t directions are indistinguishable, therefore here / = — 2(p x +p z ), implying 
the same scheme-dependence as eq. (4.5). 

Let us now consider a homogeneous system exposed to the external field, with the thermodynamic 
potential proportional to the three-volume, T log Z(Li,T,eB) = V ■ Q(eB,T). Such a homogeneous 
thermodynamic potential describes, for example, free charged particles (see, e.g., ref. [9]). On the 
one hand, the volume only appears as the proportionality factor so that in the £>-scheme, homo- 
geneity implies isotropic pressures. On the other hand, expressed as function of the magnetic flux, 
TlogZ(Li, T, <I>) = V ■ Q(<&/(L x L y ),T), which is not proportional to the volume anymore. Thus, 
in the ^-scheme, one obtains anisotropic pressures. The difference between the pressures in the two 
schemes is fixed by eq. (4.4): 



Px. 



$=Pz, pj$=p z -M.eB. (4.6) 



The two schemes give different transverse pressures, since they result from a different compression of 
the system: either pushing the magnetic field lines together with the box in which the system resides 
(<3?-scheme), or leaving the field lines unaffected (5-scheme). In the ^-scheme, one may therefore 
think of the field lines as being frozen in. This, for instance, is generally assumed to be the case for 
the magnetohydrodynamics of a perfectly conducting plasma [1, 27]. The longitudinal pressure p z is 
however independent of the scheme, since both B and <3? remain constant during the compression of 
the system in the direction of the magnetic field. 

One also has to specify the scheme when defining the pressures as the spatial diagonal components 
of the energy-momentum tensor. The tensor T Q( g can be obtained by considering the variation of the 
action with respect to the metric g 01 ^. Again, such variations can be performed keeping either B or 
<I> fixed. The conventional definition of the energy-momentum tensor (as used in, e.g., ref. [23]) corre- 
sponds to the <I>-scheme: f as a scalar quantity is not affected by general coordinate transformations, 
while fixing the two-form eB breaks general covariance. This can also be seen by considering the 
covariant derivative in a transversal direction, D x ^ y = d x ^ y ±ieBr y ^ x . Rescaling the coordinates in the 
x-y plane, (x,y) t-> (£x,£y), (d x ,d y ) \-t £,~ 1 d y ), the two terms in the covariant derivative only 
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transform in the same way (i.e., gauge co-variance is only maintained), if B transforms as B i— > £~ 2 B. 
This indeed corresponds to $ i— > Accordingly, the anisotropy, eq. (4.6), between the pressures in 
the <3?-scheme agrees with the anisotropy defined through [23], with = trT^*) (and similarly 

j(B) = t r r( B )). 

We remark that in general, besides the magnetization of the polarizable medium, logi? also 
contains the energy of the magnetic field itself, —V/T • B 2 /2. (Note that in this term the magnetic 
field appears separately and not in the product with the charge.) Therefore, in the derivative with 
respect to B, an additional term —B appears, which we have neglected above. The difference of 
the magnetic field B and the magnetization equals the external magnetic field H, which would have 
been present in the absence of the medium [48], H = B — Me. If this static magnetic energy of the 
background field is also taken into account, M ■ eB has to be replaced by (Me — B)B = —HB in 
eqs. (4.4)-(4.6). 

4.2 On the lattice 

Let us now consider the lattice setup with extents = N^a and lattice spacing a. The natural 
variable describing the magnetic field on the lattice is the flux $, which is quantized, due to the 
finiteness of the lattice volume in the directions perpendicular to the magnetic field, 

$ = q d B- N 2 s a 2 = 2nN b , N b G Z, (4.7) 

where the down quark charge qd enters. Therefore, the lattice system automatically realizes the $- 
scheme, since the flux cannot be changed continuously. In accordance with eq. (4.3), we then expect the 
lattice pressures to develop an anisotropy for nonzero B. Furthermore, we neglect the pure magnetic 
energy —B 2 /2 in logZ, since it only constitutes a constant irrelevant shift in the path integral of 
eq. (2.1), and cancels from all expectation values. This implies that the magnetic field only enters 
log J? in the renormalization group invariant combination qfB ~ eB. Altogether, the lattice pressures 
therefore satisfy 

p x -p z =Py-p z = -M-eB. (4.8) 

Due to the quantization, eq. (4.7), the magnetization cannot be determined on the lattice as 
the partial derivative of logZ with respect to the magnetic field. However, we can make use of the 
anisotropy eq. (4.8) of the <3?-scheme pressures, and compute the magnetization M from the difference 
between the transverse and longitudinal components. As shown in app. C, the anisotropics of both 
the gauge and fermionic actions enter this difference, 

- M ■ eB = -(C 9 + Cg) ■ [A(B) - A{£)\ - Q • J>(C/), (4.9) 

/ 

where Q, C, g and Q are renormalization coefficients for the anisotropic lattice parameters, which are 
of the form 1 + 0{g 2 ) in perturbation theory [49, 50], and are expected to be of 0(1) for the couplings 
used in the present study. We have seen in sec. 3.2 that A(Cf) is by about a factor of five larger 
than A(B) — A(£), and therefore dominates the right hand side of eq. (4.9). Hence, for a first estimate, 
we neglect the gluonic term and also assume (f = 1. The final step is to perform the additive 
renormalization of M ■ eB, which contains the logarithmic divergence ~ (eB) 2 loga. As discussed in 
app. A, this renormalization at T = corresponds to subtracting the total 0((eB) 2 ) contribution in 
M ■ eB. This is done by fitting the magnetization at small magnetic fields and extrapolating the result 
to zero eB: 



M r ■ eB = M ■ eB — (eB) 2 ■ lim M "^f . (4.10) 

eB->0 (eB) 2 
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The estimated renormalized magnetization is plotted in fig. 6, as a function of eB. We observe 
that after subtracting the 0((eB) 2 ) contribution, the magnetization becomes positive, indicating a 
paramagnetic QCD vacuum. The results for M r are also compared here to the hadron resonance 
gas (HRG) model prediction [9], showing a remarkable agreement for eB < 0.4 GeV 2 . Taking into 
account the independence of A(Cf) on the temperature, as discussed in subsec. 3.2, we conclude that 
the magnetization remains roughly constant up to the transition region. This is also consistent with 
the HRG description of ref. [9]. Note that in ref. [19] we presented results on the magnetic susceptibility 
of the quark spins, which indicated a diamagnetic behavior of this contribution. Here, we addressed 
the total magnetization of the QCD vacuum, which - besides the spin term - also includes the quark 
angular momentum contributions. 
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Figure 6. Estimate of the renormalized magnetization of QCD. The positivity of M r indicates that the QCD 
vacuum is paramagnetic. Shown are lattice results on five different lattice spacings (colored points) and the 
hadron resonance gas curve [9] (solid line). 

Let us summarize this section. We have shown that the pressure in the presence of an external 
magnetic field is isotropic if the U-scheme is applied, i.e., if the magnetic field eB is kept constant while 
compressing the system. Conversely, keeping the magnetic flux fixed (^-scheme) results in anisotropic 
pressures, and the difference between p x and p z turns out to be proportional to the magnetization M. 
Due to the quantization of the magnetic flux, M is not directly accessible on the lattice. However, 
exploiting the fact that the lattice implementation corresponds to the ^-scheme, we can express 
M as an expectation value of anisotropics, which can be computed directly on the lattice. The only 
complication is due to the anisotropy coefficients £ 9 , C, g and (f, which should ideally be determined non- 
perturbatively. However, a first estimate of M r can be obtained by neglecting the gluonic contribution, 
and assuming C/ = 1. The full determination of the magnetization, involving the calculation of these 
coefficients, will be performed in a forthcoming study. 



5 Conclusions 

In this paper, we have presented a comprehensive analysis of the response of gluonic observables 
of the QCD vacuum to an external (Abelian) magnetic field. We performed lattice simulations using 
dynamical staggered quarks with physical masses at several lattice spacings, temperatures and external 
field strengths. As gluons themselves do not carry electric charge, they are only affected by B indirectly, 
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through their coupling to electrically charged quarks. Still, we found several gluonic quantities that 
show a pronounced dependence on B. 

First, the gluonic contribution to the interaction measure was determined, which we observed to 
undergo magnetic catalysis at low temperatures, and inverse magnetic catalysis above the transition 
temperature. It is instructive to compare this observable to the fermionic contributions (which we 
have shown to be given by the quark condensates), exhibiting the same qualitative tendencies. This 
reflects the distinct correlation between the reaction of the gluonic and fermionic sectors to the external 
magnetic field, through the strong QCD coupling. 

Second, we found the external field to induce an anisotropy between the individual components 
of the squared field strengths, separating the chromo-electric and chromo-magnetic fields in the direc- 
tions parallel and perpendicular to B. A similar tendency is visible in the fermionic sector, where the 
anisotropy was found to be by about a factor of five larger. We have shown that the sum of these 
anisotropics (up to 0(1) factors related to lattice anisotropy renormalization) equals the QCD mag- 
netization M - an observable which is otherwise not directly accessible by standard lattice methods. 
After performing the additive renormalization of M through charge renormalization, we estimated the 
magnetization by assuming trivial anisotropy factors. We have obtained M > 0, which indicates that 
the response of the QCD vacuum is paramagnetic. We remark that paramagnetic quarks have been 
shown to imply a decreasing transition temperature based on large iV c arguments [51], in line with the 
lattice results [12, 13]. 

Third, we determined the correlator of the topological charge density for our largest magnetic 
field eB = 1.1 GeV 2 . In contrast to the above results, we found no significant difference between the 
correlators in the directions parallel and perpendicular to B. At the same time, our findings indicate 
an 0(1O -7 GeV 8 ) upward shift of the correlator at small distances caused by the magnetic field. 

To obtain qualitative insight into the response of gluonic observables to the external field, we 
examined a combined Euler-Heisenberg effective action for QED and QCD. This analysis shows that 
there is, on the one hand, no coupling of the external field to terms which contribute to the topological 
charge density and, on the other hand, that there is an anisotropic coupling to the electric and 
magnetic components of the gluon field strength tensor, with a tendency, which is in accordance with 
our numerical results. 
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A Renormalization at nonzero B 

In this appendix, we discuss the renormalization of our observables at a nonzero external magnetic field. 
In general, the logarithm of the partition function contains additive divergences in the cutoff A, which, 
on the lattice, corresponds to the inverse lattice spacing A ~ a -1 . These divergences include quartic, 
quadratic and logarithmic terms at B = [52], whereas for nonzero magnetic fields, an additional 
logarithmic divergence appears, which is canceled by electric charge renormalization [53, 54], see also 
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refs. [9, 55]. Thus, in a finite four- volume of size V4 = T/V, \ogZ has the divergence structure 

log Z = ci • V±aT + C2 • V^m 2 fcr 2 + C3 • V^rri t log a + C4 • <3? 2 log a + finite. (A.l) 

We remark that the factor accompanying the ^-dependent divergence is related to the leading coeffi- 
cient of the /3-function of the theory [9, 55], C4 • 3> 2 log a ~ • (eB) 2 log a, and the proportionality 
factor is given by the squared quark charges ^2f(qf/e) 2 , which we suppress in the following. 

Clearly, the first three divergences in eq. (A.l) cancel when computing differences A in the mag- 
netic field, see eq. (2.3). Since the ^-dependent divergence is m/-independent, one easily notices that 
all divergences cancel from the difference of the condensate mfAipfi^f. Moreover, the derivative of 
the magnetic field-dependent divergence with respect to log a at constant is already finite, making 
the interaction measure eq. (2.4) in the <3?-scheme, AI^\ also completely free of divergences. Note 
however, that the magnetic field-dependent divergence does not cancel if eB ~ <£/a 2 is kept constant, 
and thus AI^ still contains a logarithmic divergence. The magnetization of eq. (4.2) also inherits 
the logarithmic divergence. Expressing these divergences as functions of eB, we obtain 

M-eB = 2/3i-(eB) 2 logo + finite, AI^ = finite, AI {B) = -4ft • {eB) 2 log a + finite, (A.2) 

which is in accordance with eq. (4.5). For the sake of completeness, we also give the divergences of 
the energy density and the pressures. In the i?-scheme, the pressures are proportional to \ogZ. Thus, 

p'5 =p g = Pi- (eB) 2 log a + finite. (A.3) 

Using eqs. (A.2), (A.3) and (4.6), the energy density and the pressures in the $-scheme then read 

e = — /3i • (eB) 2 log a + finite, pj^ = — /?i • (eB) 2 log a + finite, p z = fi\ ■ (eB) 2 log a + finite. 

(A.4) 

These logarithmic divergences are eliminated by the simultaneous renormalization of the electric 
charge and of the bare magnetic field [9, 53-55], 

^-/3 1 (e J B) 2 loga = ^. (A.5) 

Since we are not interested in the static magnetic background energy .B 2 /2, at zero temperature this 
renormalization corresponds to subtracting the total 0((eB) 2 ) contribution from \ogZ - and, thus, 
from the observables of eqs. (A.2)— (A.4). In particular, for the magnetization this renormalization 
amounts to the procedure carried out in eq. (4.10). 

B Improvement of the interaction measure 

In this appendix, we discuss the scaling properties of the interaction measure (trace anomaly) I. 
Using its definition eq. (2.4), / can be written in terms of partial derivatives with respect to the lattice 
parameters (3, amf and However, to measure / in the <3?-scheme of app. 4, the magnetic flux is to 
be kept constant, therefore the partial derivative with respect to <I> does not contribute. Altogether, 
/ can be decomposed into gluonic and fermionic contributions 

T dlogZ v d(3 T dlogZ = d(am f ) T dlogZ 

Udloga f f ' 9 dloga'v 8(3 ' /_ dloga ' V d(am f )' [ ' ' 

such that the change in these contributions due to the magnetic field reads 

Alp = -Rf) ■ As g , AI f = -R mf ■ m/A^/ 5 (B.2) 
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where we used the definition of the gluonic action density and the quark condensate densities, eq. (2.2), 
and define the scaling factors (that are related to the lattice /3- and 7-functions) 

p d P p _ dlog(am f ) 

K P — ) K ra f — ^"j , 

a log a 1 a log a 

which have been determined non-perturbatively, along the line of constant physics (LCP) for the action 
we use in ref. [29]. We remark that the external magnetic field affects neither the lattice spacing, nor 
the LCP (see discussion in ref. [12]), therefore R@ and R mf are independent of We have shown in 
app. A that both the change in the condensate mfAtpfipf and the change in the interaction measure 
AI = A/(*) (the sum of the two terms of eq. (B.2)) are finite in the continuum limit. We remark 
moreover that both quantities are expected to approach the continuum limit as 0(a 2 ), in line with 
the scaling properties of the employed action. 

The factor R mf in front of the condensate in eq. (B.2) is related to the lattice 7-function, cf. 
eq. (19) in ref. [56], 

R mj = 1 + 7/ (a), V = 0(g 2 ) = 0(l/loga), (B.4) 

where the anomalous dimension is inherited from the mass renormalization factor Z m [57]. This 
shows that Alf (and then also AI g = AI — Alf) is also separately finite in the limit a — > 0. 
However, in contrast to the total AI, the logarithmic convergence of R mf makes Alf and AI g converge 
slowly towards the continuum limit. 

In order to improve this slow convergence, we regroup the fermionic and gluonic contributions in 
AI in the following way, 

ai = Ai g + ^2 Al f = A 4 mp + Yl A/ / mp ' ( B - 5 ) 

/ / 

A/fP = -R p ■ As g - [Rm f ~ 1) • m f A$ f 1> f , AI^ = -J^mfA^f, (B.6) 

/ / 
recovering eq. (2.5) for the improved gluonic contribution to AI. Note that the improvement correction 
is proportional to R mf — 1, which vanishes in the continuum limit, in accordance with eq. (B.4). 
Furthermore, eq. (B.6) shows that — AIj mp is given just by the renormalized quark condensate (which 
we have studied in detail in ref. [13]), implying that its convergence is improved to 0(a 2 ). Therefore, 
AI g mp is also expected to scale as 0(a 2 ). The effect of this non-perturbative improvement is shown 
in fig. 7. 



C Details on the pressure anisotropy and the magnetization 

In this appendix, we derive eq. (4.9), and show how the gluonic and fermionic anisotropies are related 
to the difference of the transverse and longitudinal pressures in the $-scheme realized on the lattice. 
Let us write down again the functional integral in Euclidean space-time, as in eq. (2.1), 

Z = f Vile- 08 * JJdet [alp + am f ] , (C.l) 
J f 

where the product runs over quark flavors /. To relate the energy density and the pressures directly 
to the partition function eq. (C.l), we need to consider anisotropic lattice spacings [49, 58]. We 
define a lattice spacing a a in direction a while a = £a Q denotes the lattice spacing in the remaining 
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Figure 7. Visualization of the effect of the improvement of the gluonic contribution to the interaction measure, 
eq. (B.6): the C(l//3) = 0(1/ log a) scaling is improved to (D(a 2 ). The red band is the continuum limit based 
on the improved data (solid points), as in the left panel of fig. 1, whereas dotted points correspond to the 
unimproved data. 

three space-time directions. The physical anisotropy - which can be determined, e.g., from the static 
potential [59] or using the Wilson flow [60] - is given by 

e = f. (c.2) 

a a 

This physical anisotropy differs from the bare anisotropics £ 9 o an d £fo, which enter the gauge and 
fermionic actions, respectively, as parameters. In the gauge sector, the total action S g - as in eq. (2.7) 
- splits into two parts, 



S, 



y 



(g -P a + ^--Pa, P* = £ £ \ Re tr P m (n), P a = ^ Retr P^(n), (C.3) 



fi<u 



with P^ u (n) being the Symanzik tree-level improved plaquettes [30], 

P,M = ~ (1 - U 2 *\n)) + | (1 - U^\n)) , (C.4) 

where uj^' n ^ denotes the product of gauge links along the boundaries of planar rectangles of size mxn 
lattice units in the \i-v plane. Note that two separate couplings (J3 and £,go(P)) need to be introduced 
in the gluonic sector to define the continuum limit along a line of fixed physical £ [49] . In the fermionic 
sector, the anisotropy shows up as an extra factor of £jq m the a direction in the Dirac operator, see, 
e.g., refs. [50, 61]. Schematically, ij) reads 

aJj> = £/o • a$ a + ^2 0,]/)^, alj)^ = - (rj^U^ - r? M L^n*^ , (C.5) 

with rjn being the staggered phases, the SU(3) links and Uu the U(l) factors, which generate the 
magnetic field B. In the end, we will set the anisotropy to £ ff o = £yo = 1 (which results in £ = 1). 
However, to calculate individual components of the energy momentum tensor, an d CgO are to be 
treated as free parameters. 
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Now, to obtain the energy density, we consider a 

1 dlogZ 



1 1 dlogZ 



V d(l/T) 

with V = (N s a) 3 being the three-volume and T 
take a as spatial, 



(C.6) 



NfN t a 4 d£ 
(Ntdt) -1 the temperature. For the pressures we 



Pc 



T dlogZ 

V a 



1 1 dlogZ 



N^Nta- 



N%a 2 a r 



(C.7) 
is the 



H u- d£ 

where L a = N s a a is the size of the lattice in the a direction and in this case V 
three- volume and T = (A^a) -1 the temperature. 

From eqs. (C.3) and (C.5), one easily obtains for the derivative of the partition function (C.l) 
with respect to the physical anisotropy £ 



dlogi? 



9{P/t 9 o) 



0( 



Pc, 



f 



0( 



+ 



d£fo 

dt: 



E 



N*N t Nr. 



(C.J 



where we used the definition eq. (2.2) of the condensates and the identities 3 







log det(alp + arrif) = Tr 



ilp c 



<9£/o alp + arrif 1 1 

The above trace is taken over space-time and color indices (note that in the staggered formulation, 
spinor indices are absent). The last term in eq. (C.8), proportional to the number of colors A^ = 3, is 
chosen such that the energy density and the pressures according to eqs. (C.6) and (C.7) are consistent 
with the definition of the interaction measure, eq. (B.l), see below. Note moreover that the parameters 
of logi? - besides /3, arrif, tgO and £/o ~ a ^ so include the magnetic flux ~ a?eB for nonzero magnetic 
fields. To obtain the pressures in the ^-scheme of sec. 4, we have to consider the partial derivatives 
at fixed <£. Therefore, the partial derivative with respect to $ does not contribute to eq. (C.8). 

Let us now consider eq. (C.8) at £ = 1. The magnetic field does not affect the line of constant 
physics for the quark masses (see the discussion in ref. [12]). Therefore, d(arrif)/d£ is independent of 
the direction a, in which the anisotropy was introduced, and the corresponding term cancels from the 
difference between p x and p z . Similarly, B has no effect on the lattice spacing [12]. Therefore, the 
derivatives of the inverse gauge couplings d(f3t g o)/dt, d(P/^ g o)/dS i and d^fo/dt are also independent 
of a. Altogether, we obtain for the difference of the pressures, 

Px~Pz = ~ \C 9 P ■ (-Px + Pz) ~ 4/5 • i-Px +P*) + (fYl ~ 
L / 

and an analogous expression for p y — p z . Above we introduced the abbreviations 



(CIO) 



Co 



1 dtftgo) 



c 9 = 



«=1 



1 dWM 



Cf 



dt 



(C.ll) 



in a manner that in perturbation theory, both gluonic coefficients Q and Qg [49] , and also the fermionic 
coefficient Q [50] are 1 + 0{g 2 ). Using eq. (C.3), and writing out explicitly the components of P^v 
according to eq. (2.8), we get, 

l ^[trfiJ(n)-tr^(n)-trfJ(n) + tr^(n)], p = x,y. (C.12) 



Po ~t~ Pz — Pa Pz 



6 ^ 

n 



3 We suppress Grassmannian integration on the right hand side of eq. (C.9), and also in eq. (C.15) below. 
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Therefore, in terms of the anisotropics eqs. (2.10) and (2.12), the pressure difference becomes 

p x - Pz = Py - Pz = ~(C 9 + 4) • 1MB) ~ A(£)} -CfYl MCf), (C13) 

/ 

as used in the body of the paper, eq. (4.9). Note that at zero temperature, the perpendicular compo- 
nents coincide, and hence A(B) — A(£) is proportional to tr£|(n) — tv B 2 (n). 

We remark that the anisotropy coefficients and the multiplicative factors in eq. (C.8) are not 
completely independent. Consider the combination I = e — p x — p y — p z , which can also be measured 
without introducing any anisotropy, as in app. B. Using eqs. (C.6)-(C.8) for the individual terms, we 
get for this combination at £ = 1, 



'-v 

where we used that 



2(C 9 ~ C 9 )P ■ (S a ) + 4 £ _ C/ £ amf 



(C.14) 



^Vflp^f = Tr a ^ = Trl - ftm/f/f/, Tr 1 = N*N t N c . (C.15) 

Note that the constant term N^NtN c /4 in eq. (C.8) is necessary to obtain a trace anomaly that does 
not contain any trivial constant shifts, so that I vanishes in the massless non-interacting case, as it 
should. The isotropic definition of the trace anomaly, eq. (B.l), must equal eq. (C.14), which leads to 
the consistency relations 

%-«,)/»—#-. i a J^pl _ (/ = _^v). (c . 16) 

a log a a£ a log a 

The sum of the gluonic anisotropy coefficients is therefore related to the lattice /3-function [49]. The 
two contributions were determined in perturbation theory for the pure gauge case, both for the Wilson 
gauge action [49] and for improved actions [62]. In the fermionic sector, the lattice 7-function shows 
up, and connects the ^-dependence of £yo and that of am,f. 

D Euler-Heisenberg effective action 

To enable a comparison of our non-perturbative results on the field strength anisotropy to perturbative 
expectations, in this appendix we compute the effective action of quarks coupled to a constant magnetic 
background field, plus constant chromo-electric and chromo-magnetic background fields to bi-quadratic 
order in the fields. This is a generalization of the celebrated Euler-Heisenberg action, in Euclidean 
space (and at zero temperature). 

The dynamics of gauge fields in the path integral is governed by the sum of the gluonic action 
eq. (2.6) and the effective fermionic action, 

S cS = -logdetC0 + m), (D.l) 

which is obtained after integrating out fermions. For simplicity, we consider here a single fermion 
flavor with mass m and charge q. For constant field strengths in QED, S e g has been computed by 
Euler, Heisenberg and Schwinger [53, 63] to all orders. For a constant magnetic field B in Euclidean 
space, it reads 4 , 

S eS (B) = -\qBm 2 e~ s coth(qBs/m 2 ), (D.2) 

87T Jq S 



4 In our conventions, the effective action corresponds to the free energy of the system, i.e. Seff = — log Z. Moreover, 
for the case of quarks, an additional factor of JV C = 3 is to be included in S e ff. 



- 19 - 



where V4 = T/V again denotes the four- volume, q is the charge of the fermion (the electron), for a 
review see [55]. This action is divergent for small s, i.e. in the UV, since coth(qBs/m 2 ) = m? /qBs + 
qBs/(3m 2 ) — (qBs) 3 /(45m e ). The leading singularity is independent of B, and is thus absent in the 
difference S e g(B) — S e Q (0). The singularity quadratic in B is taken care of by charge renormalization. 
The first non-trivial order is quartic, where, in diagrammatic language, four external photon legs 
interact via an electron loop (light-by-light scattering). The next order comes with an additional 
factor of (qB) 2 /m 4 , and is thus negligible for weak fields. 

The fourth order term for constant field strength G^ u in an arbitrary gauge group has been given 
by Novikov et al. [64], 



^ew>- 576vr 2 m 4 



{G^ v Gj_ lv ) 2 — — {G^a, G au } 2 — — [G^ajGau] 2 + — [G^u, G a /3] 2 



(D.3) 



where A is the coupling in the covariant derivative = + iXA^ and {..,..} and [..,..] denote the 
anti-commutator and the commutator, respectively. For pure QED, upon replacing \G^ U — > qF^ u , 
this yields 

^ E ' B) = -sttS t (E2 + b2)2 " 7(EB)2 ] ' (D - 4) 

reproducing the result by Euler and Heisenberg, (— E 2 + B 2 ) 2 + 7(E B) 2 , if we change from Minkowski 
to Euclidean space, by multiplying the electric field by an imaginary unit. 

In the following, we again denote SU(3) fields by calligraphic letters and U(l) fields by straight 
characters. For QCD in external magnetic fields, one has to replace XG^ — > Tnv + QB(5 fl i5 U 2 — Svifiuz) 
and a careful evaluation of eq. (D.3) to bi-quadratic order 5 yields 



0(2,2), F Vj {qBf 



3 tr B\ + tr B\ + \xS\-- tr 



(D.5) 



in terms of the field strength components defined in sec. 2.2. No topological charge term S3 appears, 
as expected from CP arguments in a purely magnetic external field. 

Thus, in perturbation theory for constant fields \qB\, \ J-^,u\ "C m 2 the chromo-electric field parallel 
to the external field has an increased action compared to the perpendicular fields, whereas the parallel 
chromo-magnetic field reduces the action. This means that parallel i^-fields are disfavored, while 
parallel S-fields are favored. This is in qualitative agreement with our non-perturbative findings that 
A{£) > and A(B) < 0. 

The remainder of this appendix is devoted to check the main formula eq. (D.5). First, let us 
revisit the Abelian theory, by removing the traces and replacing calligraphic letters by q times straight 
ones. This should be the fourth order result of Euler and Heisenberg, eq. (D.4), up to the fact that 
here we have split the i3-field in the z-direction artificially into B» + B and computed only the terms 
of 0(B 2 ). If we do the same in eq. (D.4), we obtain in this order 

- ^2 ^4 [( E ± + E l + B l + I^ll + B] 2 ) 2 - 7(E ± B ± + E\\ [B { \ + B}) 2 ] (D.6) 
Va a 4 o 2E l + 2E l + 2B ± + 2B l + 4B ii - 7E l 

4 q B 2 — !! ^ !! + (D.7) 



180vr 2 m 4 

which is what we get from eq. (D.5), too. 

5 We note that additional terms of the form qB tr(jBii.FiJ I/ ) and qB tr(£ 11 T^T^v) also appear, and contribute to S^j . 
For instance, eq. (D.8) below contains a term qB Ba = qB trBy. 
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For the second check of eq. (D.5), we note that to this order it suffices to consider the effective 
action of just a constant field Bp or just a constant perpendicular field, or just a constant field £11, 
respectively, plus the external field B. In such a situation, the color fields can be gauge transformed 
to a diagonal form, e.g. Bu = diag(/3i, B2, B3). Now the effective action splits into the sum over color 
sectors, in which one can use results for Abelian fields. Let us discuss the three cases separately. 

For the case of Bit, one has to replace qB — > B a + qB in eq. (D.2), with a = 1,2,3. To fourth 
order, this yields 

B) = ^ ^{Ba + qBW jH f e~ s (-^) ■ (B a + ^f, (D.8) 
and to bi-quadratic order, using ^ Q B 2 = trZ? 2 , 

sir\B r ,B) = V (%B?Bl\ = -J^¥ trBf, (D.9) 

eff ^ 11 ' ' 87T 2 ^V 2 / m V 45 / 607r 771 " 

which agrees with the first term in eq. (D.5). 

This result can be re-used for the perpendicular case. Electric fields can be turned into magnetic 
fields and vice versa, through the appropriate Lorentz-transformations. Here, we have a magnetic field 
(including charge) of qB in the z-direction and perpendicular electric fields S a in, say, the x-direction. 
A 'Euclidean boost' can be used to remove the electric field, leaving a mag netic field y/E* + (qB) 2 
in the z-direction. The magnitude is fixed by preserving the Lorentz invariant combination E 2 + B 2 
(and EB = 0). Thus, in the previous formula eq. (D.8) we have to replace B a + qB — > \/£ 2 + (qB) 2 , 
arriving at 

_ £ £ 0*>"^ B) - (D.10) 

which applies to perpendicular magnetic fields, too. 

For the £y case, one needs to consider a magnetic field qB and electric fields £ a (again including 
charge factors) in the z-direction. The corresponding effective action can be obtained in analogy to 
eq. (D.2), 

T/ POO J 

S a = ^£ a qB — e- s coth(s£ a /m 2 )coth(sqB/m 2 ). (D.ll) 
8vr 2 Jo s 

Upon sending £ a — > 0, this reduces to eq. (D.2), as it should. For the bi-quadratic order it yields 



8tt 2 



which confirms the remaining term of eq. (D.5). 
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